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1. Introduction 

Metastable helium (He*), has recently joined the list of atomic species for which 
Bose- Einstein condensates (BEC) have been realized jj], |j. Its major particularity 
is the 20 eV internal energy of the metastable state. Although this metastability 
leads to additional possible loss channels, it has been shown that these are not a 
serious problem. Indeed, ionizing collisions are a benefit because their low rate is 
nevertheless easily detectable. Ion detection is thus a new, "non-destructiv" and 
real-time observation tool for studies of phenomenon of BEC formation kinetics 
H [|, ||, |[ 0] ■ In this paper we will describe our progress toward rendering quantitative 
the ion signal. 

Several loss mechanisms are specific to the metastable state. First, collisions with 
the background gas lead to Penning ionization of the background gas: 



The positive ion X + produced can be easily detected and if this is the dominant ion 
production mechanism as it is for a dilute sample (for a density n < 10 10 cm -3 ), the 
corresponding flux is proportional to the number of trapped He* atoms. So for example 
we can easily measure the lifetime of a dilute, trapped sample. This linearity no longer 
holds when the density of the trapped cloud becomes high. Collisions between atoms 
in the cloud must be taken into account. The relevant ionization mechanisms involve 
both two-body processes: 



When these processes are present, the ion flux is related to the spatial integral of n 2 
and n 3 . At BEC densities, the 2- and 3-body processes dominate the background gas 
ionization, and so detecting the ion flux in this case amounts to monitoring the atomic 
density. 

In this paper, after a rapid description of our experimental setup, we present 
observations, via the ion flux, of the formation and the decay of a He* BEC. The 
observations are mainly qualitative, but we discuss some of the requirements for 
making them quantitative. We then discuss our measurements of the 2- and 3-body 
ionization rate constants both in a BEC || and in a thermal cloud. We discuss some of 
the systematic errors in these measurements and conclude with some ideas for avoiding 
these errors. 

2. Setup and experimental procedure 

Our setup has been described previously |l], |[ |9|. Briefly, we trap up to 2 x 10 8 
atoms at 1 mK in a Ioffe-Pritchard trap with a lifetime (r) of 90 s. We use a 
"cloverleaf configuration (Fig.|l|) |n| with a bias field B a = 300 rnC. The axial 
and radial oscillation frequencies in the harmonic trapping potential are typically 



= 47 ±3 Hz and v±_ = 1200 ±50 Hz respectively (w/2?r = (v^l) 1 / 3 = 408 Hz). In 



X + He* -> X + + He + e 




(1) 



as well as a three-body process: 

He* + He* + He* 



He* 2 +He*{~ ImK) 
^ He+ + He(lS) + e 



(2) 
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Figure 1. Experimental setup. The cold atoms are trapped in a cloverleaf type 
magnetic trap. A special feature of our set up is the microchannel plate detector 
(MCP) placed below the trapping region. Two grids above the MCP allow us 
either to repel positive ions and detect only the He* atoms suddenly released from 
the trap (time- of -flight measurements), or to attract and detect the positive ions 
produced in the trapped cloud (ion rate measurements). 



a typical run, forced evaporative cooling takes place for 40 s and is divided into 4 linear 
ramps. The last ramp lasts 5 seconds and the frequency decreases from 2000 kHz to 
a value between 1500 and 1000 kHz, depending on the condensed fraction wanted. A 
frequency of 1000 kHz (which is about 50 kHz above the minimum of the trapping 
potential) corresponds to the formation of a pure condensate. 

A special feature of our set up is the detection scheme, based on a 2 stage, single 
anode microchannel plate detector (MCP) placed 5 cm below the trapping region 
(Fig.Q). Two grids above the MCP allow us either to repel positive ions and detect 
only the He* atoms, or to attract and detect positive ions produced in the trapped 
cloud. To detect the ion flux, the MCP is used in counting mode || : the anode pulses 
from each ion are amplified, and processed by a counter which records the time delay 
between successive events. We can also use the MCP to record a time-of-flight signal 
(TOF) of the atoms released from the trap. Because the width of the TOF distribution 
is small (about 5 ms for a BEC) compared to the mean arrival time (100 ms), all of 
the atoms hit the detector with nearly the same final velocity of 1 m/s. The TOF 
spectra are then proportional to the spatial distribution along the vertical direction, 
integrated over the two horizontal directions. To record the TOF we use the MCP in 
analog mode to avoid saturation due to a very high instantaneous flux S . 
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3. Monitoring the evolution of a He* cloud 

To monitor the evolution of an atomic cloud, one usually releases the cloud and 
measures the TOF signal. Such a technique is destructive, and one must repeat 
the cooling sequence for each measurement. The TOF signals are thus subject to 
fluctuations in the initial number of atoms. In our case, we have a supplementary 
signal: the ion rate. We can thus minimize these fluctuations, by selecting runs 
having identical ion rates from the time between the beginning of the last rf-ramp 
until release. 

Another type of observation is possible, however. We can use the evolution 
of the value of the ion rate, which is obtained in a single run, independent of any 
initial fluctuations. When the density is close to the density for BEC formation 
(i.e. n > 10 12 cm -3 ), 2- and 3-body collisions within the cloud dominate the ion 
production. Thus the ion rate is related to the density of the cloud via the 2- and 3- 



body rate constant. Under some conditions (see Appendix A) a record of the ion rate 
followed by a TOF measurement at the end of the formation of the BEC allows one 
to monitor the evolution of all the parameters of the cloud. In such an observation, 
knowledge of the 2- and 3-body rate constants is essential. This is the aim of the 
experiments described in section |[ 

3.1. Observation of condensate formation during the evaporation ramp 

Before trying to do a quantitative experiment on BEC formation out of a non- 
equilibrium uncondensed cloud |j, Q], we can explore qualitatively what happens 
during our standard evaporation ramp. We show in Fig. || the evolution of the ion 
rate from 2 seconds before the end of the rf-ramp to 2.5 seconds after it. In addition 
we show the TOF signals corresponding to various times before the end of the ramp, 
selected using their initial ion rate. Between times t — —2 s and t — 0, the rf-frequency 
was ramped down linearly from 1.4 MHz to 1 MHz. At t = a pure condensate is 
formed. The comparison of the TOF and ion data first shows that the appearance of 
a narrow structure in the TOF spectrum corresponds, as closely as we can observe it, 
to an abrupt change in the slope of the ion signal. Thus, not only is the ion signal a 
reliable indicator of the presence of a BEC, but also a precise measure of the time of 
its appearance. 

One's first reaction in looking at the ion rate signal is to assume that the higher 
the ion signal, the larger the BEC and the smaller the thermal cloud. Fig. [2] shows 
however, that this is not quite the case: the maximum in the ion signal arrives before 
the achievement of a pure BEC. In fact, computing the value of the ion signal is 
rather complex. First, as was discussed in Refs. (ll], [l^], as well as below, the 
indistinguishability of the atoms in the BEC renders the effective 2- and 3-body 
collision rate constants smaller than in the thermal cloud by factors of 1/2! and 1/3!, 
respectively. Collisions between condensed and non-condensed atoms must also be 
taken into account Jll| and the degree of overlap between the two clouds must be 
calculated. Thus it might be conceivable to see the ion rate go down when a BEC 



is formed. We show however in Appendix A that for a fixed total number of atoms, 



the ion rate monotonically increases as a BEC becomes more and more pure. The 
observation in Fig. || is explained by the fact that, up until t = in Fig. |[ the atoms 
are being evaporatively cooled, as well as undergoing ionizing collisions and thus the 
total number of atoms must be decreasing. An explicit calculation including the atom 
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Figure 2. Single shot measurements of the ion rate versus time and the 
corresponding TOF signals. Forced evaporative cooling takes place until t = 
(only the last 2 seconds of the RF-ramp are shown : from 1400 kHz to 1000 kHz). 
The arrows indicate the time the trap was switched off to record the TOF. The 
black curves superimposed on the TOF signals are Gaussian fits to the wings of 
the TOF. See also caption of Fig. Id 



loss is given in Appendix A and agrees qualitatively with our observations. 



3. 2. Observing the decay of the condensate 

Figure || shows a series of TOF spectra taken after the end of the rf-ramp. Two 
situations are shown. In one case (grey ion curve) the rf-knife was held on at the 
frequency corresponding to the end of the ramp. In the other case (black ion curve) 
the rf-power was turned off completely at the end of the ramp. The data show that the 
condensate remains pure with the rf-knife kept on. In the absence of the rf-knife, the 
ion rate decays much faster and one sees that the sample rapidly acquires a thermal 
component. Since the total number of trapped atoms in the presence of a knife must 
be smaller than or equal to that in the absence of rf-knife, we conclude that the rapid 
decline in ion rate is due to a loss in sample density and not in the total number of 
atoms. This conclusion is confirmed by a fit to the thermal wings which reveals a 
heating as shown in Fig. 

3.3. Measuring the total number of atoms. 

An attempt to measure the total number of atoms as a function of time is shown in 
Fig. H Both the total number and the condensed number as derived from fits to the 
TOF signals of Fig. || and Fig. || are plotted. Surprisingly the total number of atoms 
appears to increase between t — and t = 1 s. There must be a systematic error, 
which we can account for by recalling that in our apparatus we only detect atoms 
which make non-adiabatic transitions to the (field insensitive) m — state during 
the turn-off of the magnetic trap JD. The fraction we observe is of order 10%. It is 
quite possible that this non-adiabatic transition does not occur with equal probability 
at every point in the trap. Thus clouds with different spatial distributions may be 
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Figure 3. Same as the previous figure except that we examine the decay of the 
ion signal after t = 0. The upper TOF curves correspond to the upper (grey) ion 
decay curve (rf-shield present). The lower TOF curves correspond to the lower 
black ion curve (without rf-shield). This shows that the rf-shield is maintaining 
a quasi-pure BEC during the decay, and that in the absence of an rf-shield the 
condensate rapidly heats up, causing the ion rate to drop even faster. 
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Figure 4. Heating of the condensate in the absence of an rf-shield. The 
temperature increases from 1.1 /iK to 2.2 fiK in 1.5 seconds. The time t = 
is the same as in Fig. ji| and |J. For each different time, 4 different TOF's have 
been acquired and fitted. 



converted to the m — state with different efficiencies. This could explain why atoms 
in the thermal cloud are observed with a higher efficiency than condensed atoms, as 
indicated in Fig. ^. 

We conclude that our measurements of the absolute number of atoms contain 
uncontrolled systematic errors of the order of a factor of two. So, even if we knew 
the ionization rate constants, we cannot use the ion rate to study condensate growth 
kinetics because we need the absolute value of the initial number of atoms and it would 
also be useful to measure the variation of the number of atoms during formation. Such 
a study will have to wait for a more reliable method of releasing the atoms from the trap 
(see conclusion). However, the measurement of the ionization rate constants is a first 
step. For a BEC, we can circumvent the systematic error on the detection efficiency 
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Figure 5. The measured number of atoms as a function of time. Crosses 
represent the total atom number, circles represent the number of atoms in the 
condensed part. The data comes from the fits of the TOF's presented in Fig. a 
and Lj and corresponds to the case where the rf knife is absent. The time scale 
indicated is the same as in Fig. |^ and |^. The increase in the total number after 
t = is spurious (see text). 



of the atoms to make a measurement of the ionization rate constants. This has been 
described in Ref. |8| and will be summarized in the following section. Afterwards, we 
will investigate the effect of this systematic error on the measurement for a thermal 
cloud. 



4. Rate constants of ionizing collisions 

The usual method of measuring the inelastic rate constants relies on fits to a non- 
exponential decay of the number of atoms. This method has some practical problems 
if the sample heats during the measurement: the density changes and complicates the 
fitting procedure. A way to avoid this heating is to apply an rf-shield, but this latter 
causes atom losses, which are not due to collisions. What is even more inconvenient 
in our case, is that what is measured in this kind of experiment is a decreasing atom 
number due to losses, which can be due to ionizing as well as non-ionizing collisions. 
We want to relate the ion rate to the density of the cloud, so what we need is the 
rate constants for ionizing 2- and 3-body collisions. We therefore use another method 
which consists in directly observing the products of the collisions, namely the number 
of ions, as a function of the density of the cloud. 

As we have seen in section ||, there is a systematic error on the measurement of 
the number of atom and so of the density of the cloud. But we will see that we can 
circumvent it in the case of the BEC. Let us then assume in a first part that we are 
able to measure accurately the number of atoms. 

We use the MCP to detect both the ions and the TOF signal. In a single run 
we record the ion rate during the last seconds of the ramp until we switch off the 
magnetic trap and record the TOF signal (to obtain the atom number N and the 
density) . The very last value of the ion rate recorded corresponds to ions produced by 
the cloud observed with the TOF signal. We repeat this sequence many times with 
different numbers of atoms in the cloud. The way to vary this number is to keep the 
atoms in the trap with an rf-shield kept on. In this way we reduce the atom number 



and keep the temperature of the cloud constant. As explained in Appendix A, the 



relation between ion rate and density is quite complex in the case of the presence of 
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collisions between atoms in the condensed part and atoms in the thermal part. We 
therefore only examine the case of a pure BEC or a pure thermal cloud. In that case 
we can write the ion rate per atom T as follows 

ion rate 1 1 „ , , 1 r , 9 . 

=r = - + -K 2 /3(n) + - K3 L(n 2 ). (3) 

where (n) = k J n 2 dr and (n 2 ) = k J n 3 dr, n being the local density. We have 
also introduced the 2-body and 3-body ionizing collision rate constants, [3 and L, 
respectively, defined according to their effect on the density loss in a thermal gas : 
ionization — — 77 — P n 2 — L n 3 . The effective lifetime r' > r is due to ionizing 
collisions with the background gas. The numerical factors come from the fact that 
although 2 or 3 atoms are lost in each type of collision, only 1 ion is produced. The 
factors K2 and K3 take into account the fact that the 2 and 3-particle local correlation 
functions are different depending on whether it is a BEC or a thermal cloud. For 
the thermal cloud K2 = «3 = I, while for a dilute BEC, one has K2 = 1/2! and 
«3 = 1/3! |0, When the sample is very dense, quantum depletion must be taken 
into account, which modifies these factors fllj ). A measurement of (3 and L would 
allow us to test experimentally the theoretical values of K2 and K3 || . 

4-1. Rate constants for a BEC 

To determine the ionizing collision rate constants (3 and L, we need an absolute 
calibration of the number of atoms in the condensate No, and the peak density no in 
order to calculate (n) and (n 2 ). As discussed above, we do not have a good calibration 
of these quantities. In the case of a BEC however, the measurement of the chemical 
potential fi obtained by a fit of the TOF signal, gives an accurate value for the product 
no a = fim/4irh 2 , a being the scattering length. With the value of ZJ we also obtain 
the product Nqo = (1/15) (H/rriui) 1 / 2 (2/i/fiIU) 5 / 2 . Experimentally we confirm that 
u oc N^ 5 where Nd is the number of detected atoms ||. This is a good indication that 
our detector is linear and that the detection efficiency for a BEC is indeed independent 
of \i. Assuming a value of the scattering length (a = 20 ran), we have therefore an 
accurate measurement of no and Nq. We have measured the rate constants (3 and 
L for a condensate ||. We obtain by a fit to equation (^) (having corrected for the 
effect of quantum depletion and the fact that the BEC also contains a small thermal 
fraction) (3 = 2.9(±2.0) x 1CT 14 cm 3 sec" 1 and L = 1.2(±0.7) x 1CT 26 cm 6 sec" 1 . 
These values agree with the theoretical estimates [[li], |l5|| . The scattering length is 
not well-known |], |j , so we have also given (3 and L for different values of a || . 

4.2. Rate constants for a thermal cloud 

To determine the rate constants for a thermal cloud we need, as before, to determine 



the atom number and density. We cannot use the same trick as in section 4.1 to avoid 
systematic errors in the detection efficiency If we want to use the above experimental 
method for a thermal cloud we must rely on a fit of the TOF to find the atom number 



and the temperature T. In |Appcndix B| , we propose a method to determine the rate 
constants which is independent of an absolute detection efficiency, but at this stage 
we will concentrate on the same technique as used for a BEC. 

As we have shown above, the detection efficiency is expected to be different for 
a thermal cloud and we can investigate the effect of this systematic error on these 
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Figure 6. Ion rate per trapped atom (T) in a thermal cloud versus average 
density. The solid curve corresponds to the value of /3 and L deduced from the 
condensate measurements. 
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Figure 7. Same as Fig. [(] but assuming a factor of 1.5 higher detection efficiency 
of the thermal cloud relative to the BEC. The data have simply been re-scaled 
along both axes; the solid curve is the same as in Fig. ra. 



measurements. We repeat the above described experiment, this time with a pure 
thermal cloud. To begin with, we assume that the detection efficiency is the same 
for a BEC and a thermal cloud. We plot the ion rate per atom as a function of 
(n) in Fig. ^. We can extrapolate the data to obtain the vertical intercept, which 
corresponds to 1/r'. For densities corresponding to the moment of formation of BEC, 
the corresponding ion rate N/t' is negligible compared to the total ionization rate, 
meaning that we are dominated by 2- and 3-body processes (see Fig. || and ||) . To 
compare with the results obtained for the BEC, we have also plotted the curve we 
would expect using the above values of (3 and L. It is clear that the data do not agree 
with this curve. Moreover, no possible pair of (3 and L taken within their error bars 
(see ||) can transform the curve so that it agrees with the data. Nor can assuming a 
different scattering length. What can make the curve agree with the data is assuming a 
different detection efficiency for atoms in the thermal cloud. If we assume for example 
that the detection efficiency is a factor of 1.5 higher for a thermal cloud relative to a 
BEC (which is consistent with Fig.||), the curve agrees better with the data as shown 
in Fig. |. 

The dispersion of the data points is quite large. This dispersion can be understood 
by examining Fig. || in which we have plotted the same data as in Fig. [| but now 
indicating the temperature corresponding to each different point on the graph. There 
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Figure 8. Same data as in Fig. |(| but with the different temperatures indicated. 
The light grey corresponds to the point with higher temperatures (the maximum 
temperature is 5.5 (J.K), the black points with lower temperatures (minimum 
temperature 1.8 fJ-K). 



is a clear systematic variation with temperature. One possible explanation is that the 
detection efficiency is temperature dependent. This agrees with the above idea that 
the efficiency depends on the spatial extent of the cloud which is indeed related to the 
temperature. We do not know the form of the detection efficiency as a function 
of temperature, but comparing these data (indicating that cold atoms are better 
detected) with the fact that a thermal cloud is better detected than a BEC, leads 
us to conclude that there exists a certain temperature giving a maximal detection 
efficiency. Therefore the correction to the detection efficiency for thermal atoms is 
not just a simple factor, but rather a function of temperature. Without knowing 
this correction, we cannot use this method to determine the collision constants for 
a thermal cloud. Still, these results are a consistency check on the rate constants 
measured using a BEC. 

5. Conclusion 

We have seen that the benefits of ion detection are twofold. First, the ion rate can 
be used to select BECs' with very similar parameters out of a sample with large 
fluctuations. Second, the ion rate itself can give information on the condensate on a 
single shot basis. Quantitatively, we still have some difficulties interpreting the data 
due to systematic errors in the detection calibration. 

One way to overcome this problem is to release the atoms from the trap by the 
mean of Raman transitions. It should be possible to transfer close to 100% of the 
atoms into the m = state. This will eliminate the temperature dependence of the 
detection efficiency and allow us to obtain more precise measurements of /3 and L, 
both for the BEC by improving the value of the scattering length and for the thermal 
cloud by making the detection efficiency temperature independent. 
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Appendix A. Predictions of the ion rate during the formation of the BEC 

The 2- and 3-body ion rates {I 2 b and /3b, respectively) in a sample containing both a 
BEC and a thermal cloud is given by |ll| : 

hb = I / dT [ n " (r) + 4n °( r ) n *^ r ) + 2n2 th(*)} (AT) 

hb = I / dT [ n " (r) + 9n o( r Hh( r ) + 18n (r)n^(r) + 6nf fc (r)] (A.2) 

where no(r) is the local density of the BEC and nt^(r) is the local density of the 

thermal cloud. Here we have taken into account the symmetrization factors, but 

neglected quantum depletion. 

Four parameters are needed to determine the densities of the two clouds : No, fi, 
and T t h- In the Thomas-Fermi approximation however, the BEC density depends 

only on /x: 

fi-U(ry 



n (r) 



0, 



(A.3) 



with U(r) the harmonic trapping potential and g = 47r ^ - the interaction strength. 
The density of the thermal cloud depends on two parameters. But, if thermodynamic 
equilibrium is reached, taking into account the interactions between the BEC and the 
thermal cloud (and neglecting the interaction energy of the thermal cloud), we can 
write: 

^(r) = -^ 53/2 (exp-^M+ 2 -°^) (A.4) 

+00 n 

where XdB is the thermal de Broglie wavelength and g 3 / 2 ( x ) — X) "^72- I n that case, 
given fj,, n t h only depends on one additional parameter. 

Appendix A.l. Comparison on the ion rates created by a thermal cloud at T — Tc 
and a pure BEC 

Before trying to calculate the ion rate for any T, which requires numerical calculation, 
let us first examine the ion rate created by a thermal cloud at T — Tc with a number 
of atoms N and that created by a pure BEC (T = 0) with a number of atoms nN 
(tj<l). 

In the case of 2-body collisions, the ratio i?2b of the ion rates created by a pure 
BEC (I BEC ) and by a thermal cloud (I th ) is related to the ratio of the peak densities. 
For 3-body collisions the ratio (i?3&) is related to the square of that ratio. Using the 
above equations we find: 

^)=C x v 2 / 5 xN- 1 / w (°) 3/5 (A.5) 
nth J ' \aj 

*2b \ a J 
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rBEC 6/5 

= %r = c 3 x ^ x r- ) (A.7) 

where a = yj -^=. The numerical factors Co — 0.78, C2 — 1.05 and C3 ~ 0.49 are 
independent of the atom considered and only assume that the cloud is trapped in a 3D 
harmonic trap. The maximum ratios are reached in the case of no loss (rj = 1). Using 
the typical values of our experiment (a ~ 20 nm, N ~ 4. x 10 5 , and uJ ~ 2tt x 408 Hz), 

we find ( ^- ) ~ 4, (^tjmax - 5 and (R 3b ) ma x - 12. 

V / max 

If the total number of atoms decreases during the formation of the BEC, these 
ratios rapidly fall. For instance, if the number of atoms decreases by a factor of 3.5 
during the last 750 ms of evaporation as shown in Fig. ^[ we would not have seen 
an increase of ionization rate but roughly the same ion rate at t = —750 ms and at 
t = ms ! This is an additional evidence of the difference of neutral atom detection 
efficiency for a thermal cloud and BEC (i.e. the total number of atoms decreased by 
less than 3.5). 



Appendix A. 2. Evolution of the ion rate between T — Tc and T = 



Using equation (Al) and (A2), we have numerically calculated the ion rates for all 
temperatures. If the cloud is at thermodynamic equilibrium all the parameters of 
the cloud are deduced from 2 parameters, for example the total number of atoms 
and the temperature. To simulate a time evolution of the ion rate we thus need a 
model for the variation of these parameters. In this appendix we will assume a linear 
evolution of the temperature between T = Tc and T ~ in 0.7 sec. This is of course 
a simplification, but given the linearity of the evaporative cooling ramp, it is a quite 
good approximation. 



In Fig. Al a) we show the evolution of the ion rates assuming a constant total 
number of atoms. The ion rate increases monotonically. We also see that the number 
of ions produced and thus also the number of lost atoms is not necessarily negligible 
compared to the total. 

We can attempt to take into account these losses in our model. In the experiments 
described in the text, the losses are not only due to the ionizing collisions but also 
to the rf-knife. In addition, losses not only lead to a decrease of the total number 
of atoms but also to a change of the temperature because these collisions change the 
condensed fraction. Thus, modelling the ion rate can be quite complicated. Here we 
wish simply to illustrate the effect of loss, so we assume that losses are only due to 



ionizing collisions, and we will neglect losses due to the rf-knife. Figure Al b) shows 
the results. The atom number decreases by only 30 % and the ion rate reaches a local 
maximum before the formation of the pure BEC, as in our experiment. Extensions 
of our model to include losses due to the rf-knife would allow one to monitor all the 
parameters of the cloud using the ion signal. 



Appendix B. Proposed measurement of rate constants independent of 
absolute neutral atom detection efficiency 

We will assume in this section that the absolute ion detection efficiency is known, and 
that 2- and 3-body losses are ionizing collisions |Q. The idea behind this method is 
that two TOF signals separated by a given time can measure the relative atom loss 
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Figure Al. Evolution of the ion rate with time. In a) the total number of atoms 
is constant. We show the different contributions to the total ion rate (continuous 
curve) of the ion rate created by 2-body collisions (short dashed curve) and 3-body 
collisions (long dashed curve.) In b) the total number of atoms (dashed curve) 
decrease with the only effect of ionizing collisions. Here the ion rate exhibits a 
maximum before the formation of a pure BEC. For both graphs, the initial number 
of atoms is 6 X 10 s and a linear evolution of the temperature between Tq and 
is imposed. The rates have been calculated with the values of rate constants 
measured in [pi. 



during this time, while the ion rate can measure the absolute atom loss. These data 
allow one to extract the rate constants without relying on an absolute calibration of the 
neutral atom detection efficiency. The method works if the neutral detection efficiency 
is unknown, but independent of temperature. Otherwise, we must also assume that 
the cloud does not heat during the measurement or that we know the variation of 
detection efficiency with temperature. 

To simplify the discussion we will neglect 3-body reactions and assume that the 
sample does not heat during the measurement. This will allow us to derive analytical 
expressions, but the results are easily generalized to include heating as well as 3-body 
reactions. We can then write the ion rate I(t) as: 

I(t) = ^ + #-"(*) a (B.l) 

T 2V eff 

with t' the life time due to ionizing collisions, N(t) the absolute atom number, V e /f 
defined by (n) — and e the ion detection efficiency. We write Nd{t) — aN(t) 
where N^it) is the detected number of atoms, and a is the neutral atom detection 
efficiency. Then 

/w - ^ + (b.2) 

OiT a ZVeff 

We can also write an equation for the atom number 

^m = -m.j. N ^ (b.3) 

dt T V e ff 

with t the total life time of the sample that we can measure independently at lowest 
density. The solution is: 

N(t) 1 



N(t ) (l + _g_j\r(t )T)e('-*o)/T - y^N(t )T 



(B.4) 
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substituting again N d (t) = aN(t) we have: 

Nd(t) = 1 fB 5) 

N d (t ) {i + J- Nd{ t Q )r)e^oVr-J- Nd {t ) T 

Thus we can measure an initial ion rate and the corresponding detected atom number 
Nd(t()) by a TOF signal, let the system evolve during a certain time and then again 
measure the ion rate and the atom number Nd(t). With the evolution of the ion rate, 



we can deduce e/ar' and — from equation (B.l), and from the evolution of the 



atom number we can deduce a y using equation (B.5). With the value of V e ff and 



e, we can obtain the value f3. We can also obtain the detection efficiency a. 

If we allow for three-body reactions, the method can still be used but (B^) is no 
longer analytical and must be integrated numerically. If the sample heats during the 
measurement, we only have to recalculate the volume V e ff for each TOF measurement. 

The reason why we have yet not been able to apply this method is as indicated 
above that the sample is heating so that the detection efficiency changes during the 
measurement. As we have not been able to measure the temperature dependence of 
a{T) the above equations cannot be solved. We hope to render the detection efficiency 
temperature independent in the near future by using Raman transitions as mentioned 
in the conclusion. 
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